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Abstract. Motivated by Harper Hamiltonians on skeletal graphs 
and their C*-geometry, we study a certain class of graph Hamilto- 
nians. These Hamiltonians can be thought of as a finite groupoid 
representation in separable Hilbert spaces. Here the groupoid is 
the path groupoid of a finite graph. Given such a setup, we con- 
sider the possible matrix versions of the Hamiltonian, which are 
indexed by the choice of a rooted spanning tree and an order of 
the vertices. The first result is that all the matrix representations 
arc linked to each other via the conjugation action of a re-gauging 
groupoid. 

We furthermore show that the symmetries of the underlying 
graph give rise to an action on the Hamiltonians of a group of 
extended symmetries. The new concept for the extension is to 
allow phase transformations on the vertices. In the commutative 
case, we prove that the extended symmetries act via a projective 
representation giving rise to isotypical decompositions and super- 
selection rules. 

We then apply these results to the PDG and honeycomb graphs 
using representation theory for projective groups and show that all 
the degeneracies in the spectra are consequences of these enhanced 
symmetries. This includes the Dirac points of the Gyroid and the 
honeycomb. 



I. Introduction 

Motivated by the synthesization of a novel material [1], we started 
an analysis of Harper Hamiltonians obtained from skeletal graphs [2]. 
Abstractly this situation is captured by graph Hamiltonians. Here the 
data is reduced to a finite graph f together with a separable Hilbert 
space for each vertex and a unitary morphism for each oriented edge, 
such that the inverse oriented edge corresponds to the inverse mor- 
phism. In this situation, as we derive, there is a Hamiltonian acting on 
the direct sum of all the Hilbert spaces and several matrix incarnations 
of it, which lie in M k (srf). Here k is the number of edges and is the 
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algebra generated by the morphisms corresponding to the loops on f 
with a fixed base point. Algebraically this setup is given by a groupoid 
representation in separable Hilbert spaces, as we explain in §2.2. 

Physically the graph f arises as the quotient graph of a given (skele- 
tal) graph T C lR n by a maximal translation group L ~ Z™. Each edge 
of the quotient graph is decorated with a partial isometry operator of 
translation in the direction of any lift of that edge to the graph V. 
The Harper Hamiltonian is constructed from these isometries; in the 
simplest version (tight-binding approximation), it is simply the sum- 
mation over them. This Hamiltonian along with the symmetries of 
the given material are the main input into the noncommutative ge- 
ometry machine, which constructs a C*-algebra that encodes relevant 
information about the system out of this data. 

In the physical situation above, the algebra srf depends on the choice 
of a magnetic field. Namely, if there is a constant 5-field given by an 
anti-symmetric form G then stf = Te-the noncommutative torus- as 
we showed in [2] . This is because in the presence of a .B-field the trans- 
lations cease to commute with the Hamiltonian (and with each other) 
and hence one enters the realm of noncommutative geometry. Here 
the translations together with the Hamiltonian generate a C*-algebra, 
which we called the Bellissard-Harper algebra, SS . This approach is 
also useful in the commutative case of vanishing magnetic field. Here 
the C*-algebra corresponds to a compact Hausdorff space X, via 
the Gel'fand-Naimark theorem. This space is a branched cover of the 
n-torus T n = S 1 x • • • x S 1 , which is the space corresponding to Tq. 
One way to view this cover is as the Eigenvalues of a family of Hamil- 
tonians indexed over the torus (see [2,3]). Namely for every point in 
t E T n there is a Hamiltonian H(t) e M fc (C), where k is the number 
of vertices of f . The fiber over a the point t in the cover n : X — > T n 
is then the set of Eigenvalues of H(t). 

One main objective is to analyze and understand the branching be- 
havior or stated otherwise the locus of degenerate Eigenvalues. The 
motivation is that in solid-state physics such degenerate Eigenvalues 
may lead to novel electronic properties, as is the case, for instance, with 
the Dirac points in graphene [4]. We have found two approaches to this 
question. The first is via singularity theory, which is contained in [3]. 
The second approach, which we take here is via group theory. Namely, 
we show that the symmetries of the underlying graph give rise to an 
action on the Hamiltonian of a group of extended symmetries. The key 
observations then are that (a) non-Abelian extended symmetry groups 
by themselves can force degeneracies via higher dimensional -i.e. > 1- 
irreducible representations and (b) any symmetries, also Abelian ones, 
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give rise to super-selection rules. The latter ones can facilitate finding 
the spectrum considerably, since the Hamiltonian Eigenspace decom- 
position has to be compatible with the isotypical decomposition of the 
representation. 

In order to find such symmetry groups, we start from the more gen- 
eral setup of graph Hamiltonians. Abstractly the setting can be under- 
stood as a groupoid representation as explained in §2.2. Here a vertex 
gets mapped to a separable Hilbert space 3^ and the whole Hilbert 
space of the theory is the direct sum over all of these. The decora- 
tion of an edge is an isometry between the Hilbert space associated 
to the vertices of this edge. By fixing a rooted spanning tree, each 
edge corresponds to a loop and we can think of the isometries as just 
being unitary elements of a C*-algebra =2/, cf. §2.2.1. In the physical 
applications this algebra will be the non-commutative torus T@. In 
the commutative case, that of zero £?-field, this is T n = Tq, that is the 
C*-algebra of continuous C valued functions on the n-torus T n . 

As worked out in [2] the fixing of a matrix version of the Hamiltonian 
from the 5-field and the decorated graph f , involves a choice of a 
rooted spanning tree r of f and an order on the vertices. The choice of 
these pieces of data is, as we derive, what ultimately gives rise to the 
enhanced symmetry groups. The argument proceeds in several steps. 

We first establish that the different matrix realizations of the Hamil- 
tonian given by choosing different rooted spanning trees and orders are 
all linked by gauge-transformations — see Theorem 3.3. The specific 
gauge transformations that arise form a groupoid. This re-gauging 
groupoid Q acts transitively on the set of all the matrix Hamiltonians 
obtained from the decorated graph. We furthermore show that the 
gauge transformations can be represented as conjugation with matri- 
ces. All the matrices here are matrices with values in srf ' . In the com- 
mutative case (§3.3), we furthermore show that these matrices give a 
projective representation of the re-gauging groupoid. Just like in or- 
dinary theory of projective representations this means that there is a 
bona fide representation of a central extension of this groupoid. In the 
case of our physical applications the extension is by T n . 

If we fix a point t G T n and evaluate the matrix Hamiltonians at 
t, the re-gauging actions give a groupoid representation in matrices 
Mfc(C). Now the stabilizer groups of a particular fixed Hamiltonian 
are the sought after enhanced symmetry groups. 

For concrete applications this leaves the problem of identifying the 
points t and the stabilizer groups or at least subgroups. To address 
the latter question, we establish that the automorphism group of f 
induces re-gaugings, by pushing forward the spanning tree and the 
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orders. That is the symmetries of the graph give rise to re-gauging 
groupoid elements, actually a sub-groupoid. Going through the con- 
struction outlined above, we can restrict to this sub-groupoid and see 
that at a fixed point of the re-gauging action we get a projective rep- 
resentation of the stabilizer subgroup. 

Furthermore, in order to identify points of t where such enhanced 
symmetry groups can occur, we show that under certain assumptions, 
that hold in all cases of our initial physical interest, the classical sym- 
metries of the graph T lift to an action on the base torus T n (Theorems 
3.5 and 3.10). At points t G T n with non-trivial stabilizer groups, we 
automatically get a projective representation of these stabilizer sub- 
groups of the automorphism groups of the underlying graph, which 
commutes with the Hamiltonian. Hence we get isotypical decomposi- 
tions, which can give us non-trivial information about the spectrum 
using the arguments (a) and (b) above. 

We wish to point out that this approach is broader than that of con- 
sidering classical symmetries of decorated graphs and in the commuta- 
tive case generalizes the extensive analysis of [5], see §3.1 for details. 

We apply all these considerations to the cases of the PDG wire net- 
works and the honeycomb lattice; see §4. Our main result here is that 
at all degenerate points, which were analytically classified in [3], there 
is an enhanced symmetry group arising from graph automorphisms in 
the above way which forces the degeneracy. 

Here the G-wire network corresponding to the double Gyroid which 
was our original motivator is the most interesting case. As shown in [3] 
there are exactly two points with triple degeneracies and two points 
with degeneracy (2, 2), that is two doubly degenerate Eigenvalues. The 
automorphism group of the graph is §4. The representation theory 
becomes very pretty in this case. There are two fixed points (0, 0, 0) and 
(ir, 7r, 7r) on T 3 under the whole §4 action. The projective representation 
is just the ordinary representation of §4 by 4x4 permutation matrices at 
the point (0,0,0) which is known to decompose into the trivial and a 3- 
dimension irreducible representation which forces the triple degeneracy. 
This result was also found by [5], where an initially different system 
was considered that results in the same spectrum. 

At the point (ir, ir, ir) things become slightly more interesting. There 
is a projective representation of §4, but we can show that this projective 
representation corresponds to an extension which is isomorphic to the 
trivial extension. Hence after applying the isomorphisms we have a 
representation of S 4 and there is a trivial and a 3-dim irrep, giving the 
second triple degeneracy. The elucidates the origin of the symmetry 
stated in [5] . Notice, that the classical symmetries of decorated graphs 
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would only yield an §3 symmetry at this point, which cannot explain 
a triple degeneracy as there are no 3-dim irreps for this group. 

Things really get interesting at the two points (±7r/2, ±7r/2, ±7r/2). 
Here the stabilizer group is A4. The projective representation gives rise 
to an extension which we show to be isomorphic to the non-trivial dou- 
ble cover 2A 4 of A 4 aka. 2T, the binary tetrahedral group or SL(2, 3). 
Using the character table we deduce that the representation decom- 
poses into two 2-dim irreps forcing the two double degeneracies. 

The identification of symmetries responsible for the degeneracies at 
(±7r/2, ±7r/2, ±7r/2) is, as far as we know, new. We wish to point out 
that one absolutely needs the double cover as A 4 itself has no 2-dim 
irreps and hence the projective extension is essential. 

We also use the fact that the diagonal of T 3 is fixed by a cyclic 
subgroup C 3 of A A in order to determine the spectrum analytically. 
Here we use the super-selection rules. 

For the D and honeycomb case, we show that the degenerate points 
which are well known in the honeycomb case and were computed for 
D in [6] are all detected by enhanced symmetries. These however yield 
Abelian representations and hence we have to use the arguments of the 
type (b), that is super-selection rules, to show that the Eigenvalues are 
degenerate over these points. 



2. General setup 

We fix a unital C*-algebra stf ' . In our applications this will be the 
noncommutative torus T@. Although eventually srf will be commuta- 
tive — in particular 1g in the applications — at the start, we do not 
assume that srf is commutative. We develop the theory in the general 
noncommutative setting as far as possible for possible further applica- 
tions. 

The basic object of study will be a C*-algebra SB called the Bellissard- 
Harper algebra generated by a symmetry C*-algebra srf and a Hamil- 
tonian H . The Hamiltonian and the algebra stf are given by a groupoid 
representation. This is in a sense inverse to [2], where the algebra srf 
given by magnetic translations along a lattice was a fundamental a pri- 
ori datum. However, under certain minimality assumptions it is fully 
recoverable. 

Ultimately we will present results about the symmetries of the Hamil- 
tonian in the commutative case. The reason that these exists can 
be traced back to a projective action of a re-gauging groupoid. Said 
groupoid already exists in the noncommutative case. 
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We will first give the algebraic setting and introduce the first men- 
tioned groupoid representation. The readers not interested in this un- 
derlying layer can skip to the section §2.2.3 in which we give the matrix 
versions of the Hamiltonian. 

Before delving into the data and axioms of the general theory, we 
will briefly review the geometric constructions of [2, 6] which are our 
main motivator. Again this section can be skipped if one wishes to 
concentrate on the more formal aspects, although it illustrates the way 
these structures arise and are interlinked in a very good manner. 

2.1. The geometric setup of wire networks. The PDG examples 
are based on the unique triply periodic CMC surfaces where the skeletal 
graph is symmetric and self-dual. Physically, in the P (primitive), 
D (diamond) and G (Gyroid) case, one starts with a "fat" or thick 
version of this surface, which one can think of as an interface. A solid- 
state realization of the "fat" Gyroid aka. double Gyroid has recently 
been synthesized on the nano-scale [1]. The structure contains three 
components, the "fat" surface or wall and two channels. Urade et 
al. [1] have also demonstrated a nanofabrication technique in which 
the channels are filled with a metal, while the silica wall can be either 
left in place or removed. This yields two wire networks, one in each 
channel. The graph we consider and call Gyroid graph is the skeletal 
graph of one of these channels. The graph Hamiltonian we construct 
algebraically below is the tight-binding Harper Hamiltonian for one 
channel of this wire network. The 2d analogue is the honeycomb lattice 
underlying graphene. Graph theoretically the quotient graph for the 
honeycomb is the 2d version of that of the D surface, but as we showed, 
the behavior, such as the existence of Dirac points [3] , is more like that 
of the Gyroid surface. 

To formalize the situation, we regard the skeletal graph of one chan- 
nel as an embedded graph r C 1". The crystal structure gives a 
maximal translational symmetry group L which is a mathematical lat- 
tice, i.e. isomorphic to Z n , s.t. f := T/L is a finite connected graph. 
The vertices of this quotient graph are the elements in the primitive 
cell. 

The Hilbert space Jt? for the theory is i 2 (V(T)), where V(r) are 
the vertices of T. This space splits 0„ e y(f) where for each vertex 
v E V(f), J4? v = e 2 (n- 1 (v)) where tt : T ->■ F is the projection. All 
the spaces J^ v are separable Hilbert spaces and hence isomorphic. Fur- 
thermore if e is a directed edge from v to w in F then it lifts uniquely 
as a vector to lR n , which we denote by the same name. Moreover for 
each such vector there is a naturally associated translation operator 
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T-+ : ffl v — y by the usual action of space translations on func- 

tions. We also allow for a constant magnetic field B = 2ttQ where 
= Yl 6-ijdxi A dxj is a constant 2-form given by the skew-symmetric 
matrix G = If 7^ then the translations become magnetic 

translations or Wannier operators : J^ v — > ^f w \ see e.g. [7]. These 
operators are still unitary and give partial isometries when regarded on 
Jtf? via projection and inclusion, which we again denote by the same 
letter. 

The Harper Hamiltonian is then defined as 
(1) H= (U^ + U^) 

eeE(T) 

where the sum runs over all edges of f and given an edge e, e and e 
are its two orientations. Notice that U*, = U<-. 

e e 

Likewise L acts by magnetic translations. If A is a vector in L then 
U-> sends each Jj? v to itself and the diagonal action gives an action on 

A 

Jf. Since L is a lattice, the representation it generates is given by 
n linearly independent unitaries. The commutation relations among 
the Wannier operators amount to the fact that the representation is 
nothing but a copy of Tq, the non-commutative n-torus, see [2] or [7]. 
This gives rise to a "diagonal" representation T@ — > B(Jj?) via the 
direct sum decomposition. Here, as usual, B(Jif) is the C*-algebra 
of bounded operators on Jtf?. This representation and H by definition 
generate the Bellissard-Harper C*-algebra 3§. 

To do calculations it is convenient, however, to have the Hamiltonian 
in a matrix form. This is achieved as follows. Fix a rooted spanning 
tree r of T. This is by definition a contractible subgraph which contains 
all the vertices of f and in addition one of the vertices is declared the 
root; denote it by vq. In order to obtain a proper matrix, we also need 
to fix an order < on all the vertices, where we insist that the root is the 
first in this order. For each vertex v of T there is a unique shortest path 
in r to v . This defines a choice of fixed isomorphism U„ v : ^ — > ,Jtff v 
by translations along the edge path, see §3.2 for details. 

Using these isomorphisms, we can conjugate the Harper Hamiltonian 
to a matrix Harper Hamiltonian H TtV(h< or H T < for short which acts on 
© we v(F) = ^0 if k = \V(r)l That is H Ti< = U*HU G M k (J%) 
where U = (B veV (T) ^vv an d the matrix entries which belong to Te 
act on Jf% again via translation. The "diagonal" action of T e on Jf? 
also pulls back to JifJ' via U and yields a faithful matrix representation 
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of Te. This algebra and H T < generate a faithful matrix representation 
of 38. 

Although different choices lead to isomorphic subalgebras of Mjt(Te), 
the concrete algebras depend on the choice of order and spanning tree. 
As we show, changing the underlying data gives rise to a re-gauging 
groupoid which ultimately gives rise to the symmetries forcing degen- 
eracies in the spectrum over certain points. 

2.2. Groupoid graph representations in separable Hilbert spaces 

The abstraction of the setting above is given by the following data. 
Given a finite graph f we define a groupoid representation of f, an 
association of a separable Hilbert space J{? v for each v G V(r) and an 
isometry : ,W V — > ,ffl w for each directed edge e from v to w. This 
data indeed determines a unique functor from the path groupoid of f 
to the category of separable Hilbert spaces. The path groupoid of F is 
category whose objects are the vertices of T and whose morphisms are 
generated by the oriented edges, where the inverse of a morphism given 
by e is the one given by e. Notice, that we are looking at the mor- 
phisms generated by the oriented edges, this means that Hom(v,w) is 
the set of paths along oriented edges from v to w modulo the relation 
that going back and forth along an edge is the identity. 

In particular, this lets one view 7Ti(r, v ) in two equivalent fashions. 
First as the topological ~K\ of the realization of a graph, and secondly 
as the group Hom(vo,vo) = Aut(vo) where Horn are the morphisms in 
the groupoid. 

We let &/ be the abstract C*-algebra 7Ti generates in B(Jff) via the 
representation. This is a bit subtle, as the concrete algebra depends on 
the choice of v , but of course any two choices give isomorphic algebras. 

Given a groupoid representation as above, set J$? = ©„ 6 y(f) and 
define H by equation (1). Fix a rooted spanning tree and an order on 
the vertices. Then again, we can identify Jf ~ j££ via a unitary U 
just as above. 

After having fixed an initial spanning tree the matrix Hamiltonian 
has a different description. To each edge e from v to w, we can associate 
wt(e) := Ul oV U-*Ul VQ G U{srf). That is, we can regard f as having a 
weight function on ordered edges with weights in U(srf). If e is an edge 
of r then with this definition wt(e) = wt(e) = 1. 

An alternative way of viewing this data is as a certain type of quiver 
representation, we will comment on this more below. 

2.2.1. Fundamental group and Euler Characteristic. For any finite graph 
f the Euler Characteristic is x(f) = \V{T)\ - \E{T)\ = 1 - h where &i 
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is the rank of Hi(T ) which is the same as the "number of loops". More 
precisely, if f is connected, 7Ti(f ) = F 6l that is the free group in 61 gen- 
erators. In the examples above b\ = n, the rank of L. One way to view 
a rooted spanning tree (r, v 0) is to think of it as fixing a base point v q 
and a set of symmetric generators/basis for 7ri(f,f ) = Hom(vo,vo). 
Topologically after contracting the spanning tree one is left with a 
wedge of S' 1 s. There are b\ of these one for each non-contracted edge. 
Each simple loop around one of the S 1 s gives a generator. Picking only 
one generator per loop gives a basis. 

Without doing the contraction the correspondence on f itself is given 
by all the (ordered) edges not contained in r. To each such ordered 
edge e from v to w, we associate the loop starting at vq going along 
r to v then traversing e tow and afterwards returning to Vq along r. 
Again picking both orientations gives a symmetric set of generators of 
the free group while picking only one orientation per edge fixes a basis. 
Any edge in the spanning tree corresponds to the unit, the constant 
loop. 

In this language, the weight function wt is a representation of ni 
lifted to the edges of the graph by the above correspondence. Thus, as 
long as the base points v stays fixed, the changes of spanning tree thus 
can be viewed as changing of a basis of 7r 1 (r,t> ). If v Q moves, say to 
v' , then, as usual, any path from v' to Vq gives an isomorphism taking 
7Ti(f , v'q) to vri(f , t> ). Both types of isomorphisms will play a role later 
in the symmetry group actions. 

2.2.2. Non-degeneracy and toric non-degeneracy. We call a groupoid 
representation non-degenerate, if the images of basis of the free group 
given by the construction above are independent unitary generators of 
srf and call it toric non-degenerate if stf is isomorphic to the noncom- 
mutative torus T@. 

Notice that if =2/ is commutative and non-degenerate, then srf ~ T bl , 
the C* algebra of the torus of dimension b\ . 

2.2.3. Hamiltonian and stf from a weighted graph. After the prepara- 
tion step above, we end up with a Hamiltonian and symmetry algebra 
srf as above. One can also choose this as the starting point, which is 
what we will do now. 

Fix a finite connected graph T, a rooted spanning tree r of T, an 
order < of the vertices of T such that the root of r is the first vertex, a 
unital C* algebra £/, and a morphism wt : {Directed edges of T} — > srf 
which satisfies 
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(1) wt(e) = wt(e)* if e and e are the two orientations of an edge 
e. 

(2) wt(e) — 1 e srf if the underlying edge e is in the spanning tree. 

In general, if wt is as above and it satisfies the first condition, we will 
call it a weight function (with values in srf) and if it satisfies both 
conditions, a weight function compatible with the spanning tree. 

Fix a separable Hilbert space J^ . By Gel'fand Naimark repre- 
sentability we realize srf C B(J^ ) and we shall use this representa- 
tion. 

We shall also postulate that srf is minimal, that is that it is the 
C*-algebra generated by the wt(e) where e runs through the directed 
edges of f . This makes the terminology of §2.2.2 applicable. Also, we 
see that this is again just a lift of a representation of 7Ti(f,t> ) to the 
edges of T using the spanning tree r. 

Given this data, let k be the number of vertices of T. We will enu- 
merate the vertices v , . . . , i>fc_i according to their order; v being the 
root. Given this data, the Hamiltonian H = H(T,t, <,w) is the k x k 
matrix whose entries are in si 

(2) H l3 = wt & 

directed edges e from Vi to Vj 

We now consider Jt? := then H = (Hij)^ is naturally a matrix 
inM fc K). 

We see that the weighted graph encodes both the Hamiltonian and 
the symmetry algebra srf . 

Remark 2.1. In the general noncommutative case, this is not quite 
enough for the whole theory, as we do not recover the action of srf on 
M' and the connection between the action of H and that of &f . Recall 
that the action of srf on ffl = J$? v * = ® veV (T) ^ s gi ven on each 
summand J4? VQ corresponding to v by pulling back the action from Jt? v . 
That is, the action is a conjugated action. 

In the commutative case this is not an issue as the representation is 
exactly the diagonal representation. 

2.2.4. Geometry in the Commutative case. If SS is commutative (and 
hence also as for instance if O = 0, then there is a geometric 
version of these algebras which can be understood as the spectra of a 
family of Hamiltonians over a base. We have the following inclusion 
of commutative C* algebras % : srf SB, by Gel'fand Naimark this 
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gives us a surjection of compact Hausdorff spaces it : X — > B where 
C(B) ~ srf and C(X) ~ 3§. The correspondence is given via charac- 
ters. Namely a character is a C*-homomorphism x '■ ^ — > C . The 
characters are by definitions the points of B. Vice-versa any point 
t G B determines a character ev t : C(S) — >■ C via evaluation. That is 
any / G C(5) is sent to f(t) G C. Given a character x on £/, we can 
lift it to a C*-morphism x : M k (srf) — >■ M fc (C) by applying it in each 
matrix entry. 

Thus any point t G B represented by the character x determines a 
Hamiltonian x(H) G M fc (C) via x- 

(3) (x(H)h = xiHij) 

Thus we get a family of Hamiltonians H(t) parameterized over the 
base. One can furthermore check, see [2], that tt is a branched cover 
over srf with 7r _1 (t) = spec(H(t)). 

2.3. Gauge group action. We stress that if as a matrix depends on 
all of this data. In the geometric situation of [2] this data was necessary 
to obtain a true matrix from the original Harper Hamiltonian. In the 
commutative case, changing the data does not change the Eigenvalues 
as they do not depend on the chosen isomorphisms. 

3. Symmetries 

In order to deal with the symmetries of the graph it is helpful to first 
fix the notation as sometimes questions become subtle. 

3.1. Classical Symmetries. A graph T is described by a set of ver- 
tices Vp and a set of edges E-p together with incidence relations d where 
for each edge e, (9(e) = {v, w} is the unordered set of the two vertices it 
is incident to. A directed edge is given by an order on this pair. Hence 
for each edge e there are two ordered edges by the orders (v, w) and by 
(w,v). We usually denote these two edges by e and e. The set of all 
oriented edges is called '. 

An isomorphism of two graphs V and V is a pair of bijections 
(4> V) 4>e) 4>v '■ Vr — > Vr> and 4> E '■ E r — >■ E T /. The compatibility is 
that the incidence conditions are preserved, i.e.: if (9(e) = {v,w} then 
d((f)E{v)) = {<f>v(v), (pv(w)}. Notice the also induces a map of ori- 
ented edges, the orientation of the edge 0_e(e) given by (<f>v(v), 0y(w)), 
if the orientation of e is (v,w). 

We will treat isomorphism classes of graphs from now on. Fixing an 
isomorphism class of a graph still allows for automorphisms. These are 



Both srf and £8 are unital 
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FIGURE 1. The graphs P, D and G and honeycomb to- 
gether with the preferred spanning tree and order. 

given as follows. Fix a representative of T then an automorphism is a 
pair of compatible maps (<pv, 4>e) <Pv '■ Vr -> Vr and (pE '■ E r — > E r . 

Example 3.1. Let us illustrate this for the graphs corresponding the 
PDG and honeycomb cases which are given in Figure 1; see §2.1 for 
details about the corresponding wire networks. We fix once and for all 
an isomorphism class of these graphs and then consider their automor- 
phisms using the representatives given in that figure.. 

For the P case there is only one vertex hence <py = id is the only 
possibility. However there is an S 3 action permuting the three loop 
edges. 

The D graph has the possibility of switching the two vertices and 
freely permuting the three edges. This gives the automorphism group 
Z/2Z x § 3 . The honeycomb similarly has automorphism group Z/2Z x 
§3- 

For the Gyroid, there is an § 4 worth of potential choices for <pv Now 
all these choices extend uniquely to the edges, since there is exactly one 
edge between each distinct pair of vertices and hence the symmetry 
group is exactly § 4 . 

3.1.1. Pushing forward Spanning Trees. Given a pair (r, r) of a graph 
and a rooted spanning tree, we define the action of isomorphisms and 
automorphisms by push-forward. That is an isomorphism between 
(r, r) and (r', r') is an isomorphism from T to V such that that 4> v 
maps the root of r to the root of r' and <fi restricted to r is an isomor- 
phism onto t' . 

If we have not already specified a spanning tree on V we can extend 
any isomorphism from (r, r) to it by push-forward. This means that 
we push-forward all the vertices and the edges of the spanning tree r 
to T': E T i := <Pe{E t ) and likewise push-forward the root. 

In particular, Aut(T) acts on the set of spanning trees of a fixed 
graph T. This action is not transitive in general and may have fixed 
points. 
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Example 3.2. In the cases of PDG and the honeycomb, it is a tran- 
sitive action. 

For the G graph the action is not fixed point free, there is an § 3 
subgroup fixing a given spanning tree. 

For the P graph the action is fixed point free, while for the D and the 
honeycomb although the action is transitive, there are again stabilizers. 
For the honeycomb the group fixing a spanning tree is the §2 — Z/2Z 
interchanging the two other edges, with both vertices fixed, while in 
the D case it is an §> 3 action interchanging the non spanning tree edges. 

In the triangle graph, with one edge doubled, we can see that there 
are two orbits. One where the two non-doubled edges give the spanning 
tree, there is an §2 stabilizer interchanging the non-root vertices, one 
where one of the doubled edges is in the spanning tree, there is no 
stabilizer in this orbit. 

3.1.2. Orders, Weight Functions and Isomorphisms. Also if there is an 
order on all the vertices then the isomorphisms are asked to be com- 
patible with this order and auto- and isomorphisms can be extended 
by pushing forward the order. 

A final piece of data, which we utilize is a C* algebra valued map 
wt : E^ r — > such that wt(e) = wt(e)*. In our case of interest 
srf — Tq, and even Tq later, but for now we will keep the abstract 
setting. 

Notice that a given weight function might not be compatible with 
the new, pushed-forward spanning tree, since it does not necessarily 
satisfy the condition (2). 

The idea of re-gauging is to fix condition (2) by passing to an equiv- 
alent, re-gauged weight function. 

One way to see the failure to preserve the compatibility condition is 
view the change of spanning tree as a change of basis of i\\ or rather 
an isomorphism of 7TiS via the path-groupoid. This point of view is at 
the basis of the proof of Theorem 3.5 below. 

3.1.3. Classical vs. extended symmetries. There is another natural choice 
of iso- or automorphism for graphs with weight functions. Here one 
would postulate that the weight functions are compatible. Weight func- 
tions naturally pull back via 4>Ei that is (f)* E (w)(e) = wt(0^( e ))- Using 
that (pE is an isomorphism, one can push-forward by pulling back along 

<t>E- 

One could call these symmetries classical symmetries of the weighted 
graph. These are the kinds of symmetries that were for instance con- 
sidered by [5]. 
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We will consider symmetries of the underlying graph not of the 
weighted graph. The weights are fixed by re-gauging. One way to 
phrase this is that is we utilize an extended symmetry group which 
allows for phase factors at the vertices. The details are given below. 

3.2. Gauging. One of the key observations of [2] (see §2.1) is that 
there are several equivalent matrix versions of the Harper Hamiltonian. 
This applies also to the case of matrix Hamiltonians from weighted 
graphs. 

Starting from the groupoid representation we obtain H as an endo- 
morphism of = ©„ 6 y f where the entries are partial isometries 
U vw : Jtfu — > J%,, with U wv = U~*. By fixing a rooted spanning tree 
r, one can derive matrix versions, which act on ®^, , where v is 
the chosen root. There is one more caveat: in order to have a true 
matrix, we also need to fix an order on all the vertices, which we fix 
to always start with v . Let Vi be the i-the vertex in this enumeration. 
Then we obtain a matrix H T < by using the isomorphisms U£. v which 
are defined as follows. Let vq, w±, . . . , Wk, Vi be the sequence of vertices 
along the unique shortest path from vq to V{ in r, then 

(4) U V . VQ — Uv i w k U WkWk _ 1 . . . U W2W1 U W1V0 

and Ul = (Uy )*. Given the choice of (r, <), we get the correspond- 
ing matrix Hamiltonian as 

(5) 



r >< ■ \±s v ° \jy v i ~ ~ \jy Vi KIP v o 



Of course, all the J$? Tt< are all equivalent, although not canonically. 
For the equivalence one has to choose a path from one root to the other. 
We will exploit this fact extensively below. 

In the language we have set up, we get a weight function wt(e) = 
Uv 0V UvwUw V0 and H is then the graph Hamiltonian of the representation 
which associates to each vertex the Hilbert space Jt? Vo and to each 
directed edge e the endomorphism wt(e). 

We will now consider the relationship between the different matrices 
H T < and i? T ',<'- Notice that although they are always matrices in 
Mk(^), they represent H via different isomorphisms. More precisely, 
the first and last arrow in the composition will differ. 

There are basically three situations: First, r = r' as rooted trees 
and only the order changes. Secondly, r = t' but v ^ v' and thirdly, 
the trees just do not coincide, i.e. at least one edge is different. In 
the first case the isomorphisms are simply changed by permutations 
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of the factors J^, . This means that the difference between the two 
matrix Hamiltonians is simply conjugation by M a where this is the 
standard permutation matrix. Since vq is fixed to be the first element, 
the permutation actually lies in the subgroup Sfe_i C fixing the first 
element. Sk-i then acts simply transitively on the orders. 

In the second and third case, the situation is more complicated. Of 
course the second type of change is related to an action of but things 
are not that simple, since there is a change of the space that H T i t< r acts 
on. If the tree moves, then we have to also change the weight functions 
to make them compatible with the new spanning tree. 

Taking the point of view of 7Ti the first type of transformation is 
just a permutation of the basis. Keeping the base point, viz. the root 
fixed, we get a change of basis, but when we move the base point, we 
get an isomorphic group. In doing this, we effectively use the path 
groupoid. In pure matrix language, we will show in the next section 
that the matrices change by a re-gauging, which we can realize by 
conjugation. These re-gaugings are governed by a groupoid, the re- 
gauging groupoid. 



3.2.1. Gauging in groupoid representations. As discussed in 2.2 the 
Harper Hamiltonian before fixing a spanning tree can be thought of 
as a certain type of groupoid representation. 

For such representation, we can re-gauge it by acting with automor- 
phisms of the J^ v , that is the group x ver Aut(Jtff v )) . Here J^ v is the 
Hilbert space associated to v and U-> is the morphism Jif v — > J^ w given 

by the representation where e runs from v to w. Picking an element 
in this group is the same as the assignment v (p{v) G Aut{,W v ) 
this means that U-+ gets re-gauged to ). This gives an 

equivalent representation. 

In our situation, since we have Hilbert spaces, we can look at unitary 
equivalences and restrict the automorphisms to be unitary. Also, choos- 
ing an identification of all the isomorphic separable Hilbert spaces J4? v 
with some fixed J^ VQ we can take the re-gaugings to live in the unitary 
operators on J^ VQ . 



3.2.2. Spanning tree re-gauging. 

Theorem 3.3. Given two ordered rooted spanning trees (r, <) and 
(r', <') there is a matrix M G M k (a?) with MM* = M*M = id such 
that MH T< M* = Ul,H T , <t UL n 

vqv T v vo 
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Proof. Consider the commutative diagram : 
(6) 



i=l 



vq 



v v n 



tPj <7Gtf. = Jo : ' = tPj JVv> 



V V Q 



We see that if %' = a(i) and f = a (J) so that v' v = vf 
i H r',<')i'j' = U viv'., H v'. l v'.U^y o 

= U v' v ( U v v' U v' v'., U v iV0 ) ( U v Vi H ViVj Uv jVQ ) {Ul QV . y Q ) Ul QV , o 

With fa, = U^ v .U£ nV ,U;, oVo e T@. So that if $ = ^ay(^) and M a is 
the permutation matrix of a then M = $*M CT -i and M*M = id. □ 

We choose to place M on the left of the Hamiltonian so as to get a 
left action. 

Proposition 3.4. The matrix $ above can be expressed as a re-gauging, 
that is there is an element in the gauge group whose matrix represen- 
tation is $. 

Proof. Unraveling the definition given in equation (4), we can express 
the matrix <£> by the following iterative procedure. We start at the root 
of t' and choose <f>(v r ) = id. Assume we have already assigned weights 
to all vertices at distance % from v' and let to be a vertex at distance 
i + 1. Then there a unique v at distance i which is connected to W 
along a unique directed edge e fof the spanning tree r'. Set <f)(w) = 
wt(e)(j)(v) G U(srf). If a e is the unique permutation that moves 
the order < to the order <', then if M a is the standard permutation 
matrix and let $ = diag v '. eT > (^(v'^) , then M = $*M CT -i. □ 

Of course the form of M depends on the initial choice of <j>(v o) = id, 
which amounts to using the iso , to pull-back the matrix. Any 
other choice of iso will differ by an element of srf which is then the 
value of <f> Vo . This plays a crucial role later. 

3.2.3. Commutative case. Reduced gauge group. In the commutative 
case, we can fix a character \ '■ ^ ~ ^ C and then under \ & U matrices 
become U(l) valued and all the Hilbert spaces 3^ v become identified 
with C . In this case, we can identify the gauge group action with an 
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action of £7(1) xVr on [/(l)-valued weight functions, using A = x ° <$>• 
For every oriented edge e from v to w the re-gauged weights are 

wt'(x&) ■= \(v) X (wt(e))\(w) 

Notice that we have taken the indexed un-ordered product. If we 
fix an order of the vertices, then the group acts on the vertices 
as well and the full gauge group which acts on the Hamiltonians by 
conjugation is the wreath product G — U(l) I 

We see that the constant functions A act trivially and hence to get 
a more effective action we can quotient by the diagonal U(l) action 
and consider the reduced gauge group G := G/U(l), where U(l) is 
diagonally embedded in U(l) k and E>k acts trivially. 

Abstractly U(l) k /U(\) ~ £/(i) fe_1 , to make this explicit, we can 
choose a section of G — > G. Our choice <f>(vo) = 1 is just such a choice 
of a section. The action of on the remaining k — 1 factors is then 
more involved, however. It is still a semi-direct product, but not a 
wreath product any more. This has practical relevance in the Gyroid 
case. 

The proof of the theorem above then boils down to the fact that a 
rooted spanning tree uniquely fixes a unique gauge transformation as 
follows. We let X(root) = 1 by the global gauge U(l). Now the weight 
on each vertex of the tree is fixed iteratively by the condition that 
Xwt(e) = 1. The whole set of weights then gives a diagonal unitary 
matrix and taking the product with the appropriate permutation, we 
obtain the matrix M. 

3.2.4. Gauge groupoid. Let Ham(F, &/) be the set of all matrix Hamil- 
tonian H T < with all possible rooted spanning trees with compatible 
weight functions wt : Ef — > srf . Let U{srf) be the unitary elements 
of srf . Then the full gauge group G = U(£/) k I S k acts on the orders 
and the weight function and on the Hamiltonians by conjugation and 
permutation just as above. 

Given a concrete Harper Hamiltonian, H , its set of possible matrix 
versions Ham C Ham given by the different choices on (r, <) are 
contained in the orbit of one such particular choice. In this situation, we 
want to restrict to a groupoid. Namely for each point H T < e Ham®, we 
want to only consider those elements in G which send H T < to another 
element in Ham . By what we have discussed so far, for each two 
elements in this orbit, there is an element in G which takes one element 
to the other and we did give a particular choice for it. These elements 
constitute a groupoid or rather give a representation of a groupoid 
Q. Abstractly, Q has objects pairs (r, <) of a rooted spanning tree 
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and an order of T and a unique isomorphism between any two such 
pairs. If the two pairs coincide, the isomorphism is the identity map. 
The representation is given by mapping (r, <) to H T < and the unique 
isomorphism to the morphism of Theorem 3.3. 

3.2.5. Matrix version and cocycle. In order to calculate it is prefer- 
able to work with a matrix representation of the groupoid action. The 
problem is that although the groupoid associates a matrix to each re- 
gauging, these matrices all act in different spaces. To make everything 
coherent, we have to fix once and for all a base vertex and isomor- 
phisms U VVsf : J££„ — > ,ffl v . We do this by using some initial choice of 
spanning tree (r*,^, <*) and set U VVtt := U^. 

Now, if we are given three choices if T0;<0 , H T > t< r and i? T // ;< //, let g\ = 
MU^ VQ be the morphism of Theorem 3.3 for moving from H T < to H T / t< i 
by conjugation, and likewise define g 2 , M 2 for the move from ^ to 
H T /' t< ti and g 3 , M 3 the morphism and matrix associated to the move 
from H T < to i/ r » 5< //. 

Then as morphisms g\g 2 = g 3 as one can check by direct calcula- 
tion. If we naively multiply the matrices Mi and M 2 and think of the 
product as an endomorphism of 0^ , then we are tacitly using a 
pull-back of the endomorphism given M\ on 0^ ,^f VQ . For M 2 because 
of the definition of M, we first pull-back from 0^ to . ^ us- 
ing r and then further to 0^ JJ,, that is the matrix product M 1 M 2 
represents the morphisms M o M 2 where M = U^* v 'MiUy* Vit , M 2 = 

U Z:vo U v v' M 2 U v' v U Z*ov*- Now > we S et the matrix equation M x oM 2 = 

CoM 3 with C = Ul ,Uj',,,Ul,,„ . Now C need not be the identity 

° vov v v v vq j 

in general, but it is some element in U (<$/). More precisely, we have 
U T ,U T , „U T „ G srf is a diagonal/scalar matrix. In the case of a wire 

v v' v' v' ( " vfivo to ' 

network we have that srf = T@ and the matrix above is just of the form 
all^ with ee L. 

e 

Summing up, we consider the map p : Q — > GL^g/) given by as- 
sociating to the groupoid element g moving from i? T ',<' and i? r ",<" 
the matrix M constructed in Theorem 3.3. This map is not in gen- 
eral a groupoid representation by itself. As we just calculated for two 
composable matrices gi,g 2 let c(gi,g 2 ) be the matrix representing the 
morphism C then 



(7) p(g l )p(g 2 ) = c(g l ,g 2 )p(g l g 2 ) 
Using the definitions one calculates that 

(8) c(#i, 9s) = c(g u g 2 g 3 )c(g 2 , g 3 ) 
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the value being 

c (fi'2, 92, 93) which in the analogous notation is given by U^ QV ,U^ M^ v „,U^„ t 
So technically, we have a non-commutative groupoid 2-cocycle with val- 
ues in Tq which is given by c(g, h) := MgM^M^ which is defined for 
any two composable re-gaugings (see [8] for a definition). 

One good way, albeit a bit abstract, to view the matrix multiplica- 
tion, is that we are dealing not with a representation of the groupoid per 
se, but rather, with representations of 7Ti of the re-gauging groupoid 
and our initial choice of (r*,i>*, <*) is the choice of a base point. The 
correspondence is analogous to the one discussed in §2.2.1. 

3.2.6. Re-gaugings Induced by Graph Symmetries. If we have a sym- 
metry, aka. automorphism, of the graph f, then given a fixed choice 
(t, <), we can push forward both these pieces of data to (V, <') with 
0. This means that for every H T>< any automorphism gives rise to a 
re-gauging of H T < to i? r ',<'- 

3.2.7. Lifts to Automorphisms. One interesting question for any given 
re-gauging is if there are automorphisms i[> of =2/ such that 



where, again, -0 is -0 applied to the entries. This is the type of enhanced, 
extended symmetry we will use in the commutative case. 

One way such a symmetry can arise is by a re-gauging induced by 
an automorphism is T. A stricter requirement that is easier to handle 
is that not only the matrix coefficients of the Hamiltonian transform 
into each other, but rather already the weight functions. This avoids 
dealing with sums of weights. We say a re-gauging induced by an 
automorphism of f is weight liftable by an automorphism ip of srf if 
ip(wt(e)) = wt'(<f>(e)), where wt' is the re-gauged weight function for 
the pushed forward spanning tree. 

Theorem 3.5. Given an automorphism of T there is at most one 
weight lift by an automorphism ip of the re-gauging induced by 0. On 
the generators wt(e), e not a spanning tree edge, the putative map 
is fixed by the condition if)(wt(e)) := wt?(<j)(e)), where wt' is the re- 
gauged weight function. 

Furthermore, the ijj(wt(e)) again generate srf and hence whether^ 
indeed defines an automorphism only needs to be checked on the gener- 
ators wt(e). 

Lastly, ip is induced by a base change ofir\(T). 



(9) 
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Proof. Let wt' be the re-gauged weights after moving from (r, <) to 
(t', <'). If an automorphism ip of srf that lifts exists, then it satisfies 
wt'(4>(e)) = tp(wt(e)). After fixing an orientation for each edge the 
wt(e) generate, we see that the morphism is already fixed, since by 
assumptions the wt(e) generate srf . 

In order to show that the ip(wt(e)) are generators, we will prove the 
last statement first. As discussed in §2.2.1 wt gives a representation p 
of 7Ti(r, t> ) and wt' gives a representation p' of 7Ti(F, v' ) if vq is the root 
of r and v' = (f)(v ) is the root of r', the pushed forward spanning tree. 
In t there is a canonical shortest path p T VQV > o from v' to v . Conjugating 

by this path gives an isomorphism P : 7ii(T,v' ) — > 7ri(T,t>o). This is 
in essence the definition of the path groupoid of f . Let / r (e) be the 
loop associated to e by using r as a spanning tree, see §2.2.1, then 
wt(e) = p(/ r (e)). It follows from the definition of the re-gauging that 

V(p(r(e))) = Mwt(e)) = wt\<t>&) = p'(r'(0(e))) = p{P{V' {<t>{ e)))) 

so that ip is induced by the chance of basis l T (e) — > P(/ r '(0(e))) in 
7Ti(r, v ). From this it follows that the ^(^(e)) generate. □ 

Corollary 3.6. // the groupoid representation is non-degenerate, so 
that stf is generated by the wt(e) and each non-spanning-tree edge 
gives a linearly independent generator, then the morphism ^ above is 
well-defined as a linear morphism. 

If there are no relations among the generators, e.g. in the case stf = 
T n the commutative algebra of the torus, then every automorphism is 
weight liftable, i.e. \& from above is well-defined as an algebra homo- 
morphism. □ 

We will use the corollary in §3.3 to define the enhanced symmetry 
groups in the commutative toric non-degenerate case. 

Remark 3.7. These types of symmetries might also help to explain 
the somewhat mysterious approximate symmetry between the noncom- 
mutative and the commutative case found in [6]. Here the symmetry is 
between two loci in B. Where B in the commutative case is the base for 
the family of Hamiltonians as discussed above and in the noncommu- 
tative case is the space parameterizing the background magnetic field. 
The two loci are the locus of degenerate Eigenvalues on the commuta- 
tive side and the locus of values of B where SS is not the full matrix 
algebra. From the examples PDG and honeycomb, these two loci have 
exactly the same top dimension and there are further characteristic 
features which they have in common. These considerations would lead 
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us too far astray in the present context, but we plan to return to them 
in a future paper. 

3.3. Enhanced symmetries in the commutative case. We will 

concentrate on the commutative case in the following. One physical 
feature that makes the non-commutative theory more complicated is 
that conjugating H with elements from srf usually does not leave it 
invariant. This is of course the starting point for considering the C*- 
algebra SB which contains all these conjugates. 

3.3.1. Extension. In the commutative case, U{sf) is a commutative 
group and the 2-cocycle defines a central extension Q of Q by U{srf). 
We can consider the action of this central extension, since the action 
of U(srf) commutes with the Hamiltonians, permutations and the re- 
gaugings in this case 

(10) 1 ->■ V{sf) -» 1 

Then p does give a groupoid representation of Q. 

Remark 3.8. There is a nice geometric interpretation of this in the 
case of wire networks. Here the group U(srf) corresponds to trans- 
lations along the lattice L. One can identify the vertices of f with 
the elements in chosen primitive cell and likewise one can arrange the 
spanning tree edges to be inside this cell. When we are re-gauging, we 
move the base point along the spanning tree edges. After doing this 
several times the new root can lie outside the original primitive cell. 
The co-cycle then measures the displacement of the new cell relative 
to the old cell in terms of an element A G L, more precisely it is just 
U x . 

3.3.2. Enhanced symmetry Group. In order to find degeneracies in the 
spectrum, we use the characters and then look for fixed points under 
the induced groupoid action. Using the language of §2.2.4, given a 
point x £ &4 we get a ma P X '■ Ham — » M k (C). There is then an 
induced action of the groupoid on x(_£/am ), by pushing forward with 
this character. It can now happen that x{H T ,<) — x(H T ',<'), that is 
H T ,<(t) = i? r ',<'(i), for the point t G B corresponding to x- 

For each element H(t) G x(Ham Q ), we get its stabilizer group St(H(t)) 
under the induced groupoid action. This is the image of the transitive 
action of the groupoid on the fiber of x over H{t). We can identify 
St(H(t)) with the image of that subgroupoid. If this group is not triv- 
ial, which means that the fiber is not just a point, we call this group the 
enhanced symmetry group of H(t). It is realized by re-gaugings, that is 
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conjugation by specific matrices which form a projective representation 
of the stabilizer group as we presently discuss. 

3.3.3. Super-selection rules, Projective Representation and Degenera- 
cies. If St(H) 7^ 1 then this means that the set of all matrices xp(g) 
for g G St(H(t)), where we identified g with its defining element in Q, 
all commute with the Hamiltonian H(t) and hence each one and all of 
them together give super-selection rules. This of course is already a 
great help in finding the spectrum. 

Since p is only a groupoid representation of Q, we get that xp is 
a representation of a an extension of St(H(t)). Since we are in the 
commutative case, this extension is central and gives rise to a projective 
representation of St(H(t)). 

(11) 1 -> 17(1) -> St(H(t)) -> St(H(t)) ->• 1 

In order to apply the general arguments of representation theory, we 
will be interested in the class of this extension. These extensions are 
classified up to isomorphism by H 2 (St(H(t)), U(l)) [9, 10]. 

We give a brief definition of this cohomology group, as it is impor- 
tant for our calculations (see e.g. [11]). Let G be a group and A be an 
Abelian group, which we also write multiplicatively 2 . Set C l (G,A) : = 
Map(G xl , U(l)) these are the i-th cochains. There is a general differ- 
ential d : C % — >■ C l+l with d 2 = 0. We will need the formulas for it on 1- 
and 2- cochains. If A G C 1 (G, A) then d\(g, h) = X(g)X(h)X(gh)~ 1 and 
if c G C 2 (G, A) then dc(g, h, k) = c(h, k)c~ l (gh, k)c(g, hk)c(g, h)' 1 . Set 
Z 2 {G,A) := ker{d : C 2 (G,A) ->• C 3 {G,A)) and B 2 {G,A) := Im{d : 
C 1 (G, A) — > C 2 (G, A)). Notice that an element c is Z 2 precisely means 
that c satisfies the cocycle condition (8). Now B 2 (G,A) C Z 2 (G,A) 
and H 2 {G,A) := Z 2 (G, A)/B 2 (G, A). 

What this means is that we can move to an isomorphic extension by 
using a rescaling A G C 1 (S't(iJ(t)), U (1)). Another interesting concrete 
question is if a given homology class [c] can be represented by a cocycle 
in a subgroup of A. This is especially interesting if the subgroup is 
finite. In our concrete calculations for the Gyroid, we will use for 
instance Z/2Z and this will lead us to consider double covers. 

In general if we identify that the projective action of St(H(t)) is iso- 
morphic to an action of a finite group extension St(H(t)), then we can 
use this representation to decompose C k into its isotypical decomposi- 
tions with respect to this group action. If the group is non-Abelian, 
then there is a chance that some of the irreducible representations in 



2 We consider A to be a G module with the trivial action. 
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the decomposition are higher dimensional, which implies degeneracies 
of the order of these dimensions. Again this is present for the Gyroid. 

3.3.4. Geometric lift of the groupoid action. In order to understand the 
(projective) group action, geometrically in the commutative case, one 
lifts the action on the Hamiltonians to an action on the underlying 
geometric space. We will now for concreteness fix srf — T n , that is 
the groupoid representation is commutative, toric non-degenerate, as 
is the case in all crystal examples we consider: PDG, Bravais and 
Honeycomb. 

Finding lifts then means that one considers the commutative diagram 
(12) T n 

¥ 
rpn 

where the dotted morphism is the lift to be constructed and H T< is 
the map t h- > H Tj< (t) := x(t) if \ is the character corresponding to t. 

The existence of these lifts is not guaranteed in general, and indeed 
there are examples of re-gaugings that cannot be lifted. A non-liftable 
example can e.g. be produced from the cube graph obtained from the 
Gyroid graph by quotienting out by the the simple cubic lattice, see 
[2]. We will show that all lifts stemming from automorphisms of the 
underlying graph do lift. 

Looking at the diagram (12), one consequence of this action is that it 
lets us pinpoint Hamiltonians with enhanced symmetry group. Using 
Corollary 3.6 and translating it to the geometric side, we obtain 

Proposition 3.9. Let (f , wt) be a toric non-degenerate weighted graph. 
In the commutative case the automorphism group ofT lifts via the gaug- 
ing action to an automorphism group ofT n . That is we get a morphism 
Aut(t) ->■ Aut(T n ). 

If a point t G T n is a fixed point of a lift of an element g G Q , 
then, the re-gauging is an enhanced symmetry for the corresponding 
Hamiltonian, that is x(p(g)) commutes with the Hamiltonian H Tt< (t). 

□ 

Summarizing these results: 

Theorem 3.10. If (T,wt) is commutative and toric non-degenerate, 
then a stabilizer sub-group G t of t G T n under the induced action of 
Aut(T) on T n leads to an enhanced symmetry group St(H(t)) for the 
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Hamiltonian H{t). This group also has a projective representation via 
the matrices x{p{g))- □ 

We can exploit the representation theory of this group to get infor- 
mation about degeneracies. 

Example 3.11. For commutative toric non-degenerate groupoid rep- 
resentations of symmetric graphs the re-gaugings by re-orderings are 
always representable via an automorphism of the graph. If a e Sfc-i 
permutes the vertices of the spanning tree leaving the root fixed, then 
the re-gauging lifts as the reordering of the generators and possibly 
taking * of them. The matrices are just the usual permutation matri- 
ces of C Sfc acting on the last k — 1 copies of C in C k . 

Remark 3.12. In the commutative case, the representation of 7Ti(f , v ) 
factors through its Abelianization Hi(T). 

4. Calculations and Results for Wire Networks 

In this paragraph, we perform the calculation for the general theory 
we developed for the PDG and honeycomb wire networks as described 
in 2.1. In all these situations Theorem 3.10 applies. The upshot of the 
following calculations together with the analysis of [3] is: 

Theorem 4.1. In all the examples PDG and honeycomb, all the fixed 
points come from fixed points of the semi-classical action. Moreover the 
fixed points, stabilizer groups, their extensions and the decomposition 
into irreps for the case of the Gyroid are given in Table 2. 

4.1. Gyroid. The graph f in the Gyroid case is the full square. It has 
symmetry group §4. With the Gyroid weights, the graph is faithful 
and hence the action can be lifted to an action on the torus. It acts 
transitively on all ordered spanning trees. Such a spanning tree is 
fixed by specifying a root and the order. The subgroup of §3 acts 
transitively on all orders. The matrices of this subgroup action are 
just the permutation matrices acting on the last 3 copies of C in C 4 
and the lift of the S 3 action on the generators A, B, C of T 3 is given by 
the permutation action. 

We fix an initial rooted spanning tree and order as in Figure (3). 

4.1.1. Action on T n . The action of § 4 on T 3 is fixed once we know the 
action of the generators (12), (23) and (34). 

The action of (23) is graphically calculated in Figure 2, from which 
one reads off *((23))(A, B, C) = (A*,C*,B*). Here (A, B, C) is the 
notation for the initially chosen basis of T 3 . 
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Figure 2. Calculation of the action of (23) on T 3 . The 
original graph, the pushed forward order and the move 
into the old position to read off the morphism 



In the graphical calculation, we first write down the graph together 
with the initial spanning tree and order. We then push-forward the 
spanning tree and the order. For this we keep the vertices and edges as 
well as the weights fixed. We then (if necessary) give the re-gauging 
parameters by writing them next to the respective vertices and (if nec- 
essary) perform the re-gauging. Finally we move the vertices and edges, 
so that they coincide with their pre-images to read off the morphism 
on the generators given by Theorems 3.5 and 3.10. 

A similar calculation shows that *((34))(A, B, C) = (B*,A*,C*). A 
consequence is that the cycle (234) = (23) (34) acts as *((234))(A, B, C) = 
$((23))(B*,A*,C*) = (B,C,A) and is the cyclic permutation. 

The action of (12) is more complicated as the root is moved. For this 
we calculate graphically, see Figure 3, and read off \& as: (A, B, C) >->■ 
(A*,B*,ACB). 

This allows us to compute fixed points and stabilizer groups. We 
will first concentrate on non-Abelian stabilizer groups. There are only 
two fixed points under the full S> 4 action and these are (1,1,1) and 
(— 1, —1, —1). The group At, the subgroup of all even permutations, is 
the stabilizer group of the two points (i, i, i) and (—i, —i, —i). One can 
readily check that these are the only non-Abelian stabilizer groups. 
The other possibility would be S3, but a short calculation shows that 
anything that is stabilized by any §3 subgroup is stabilized by all of §4. 



4.1.2. Representations. We collect together the matrices M needed for 
further calculation. Again, we fix our initial ordered rooted spanning 
tree as before. 

Using short hand notation, the matrix for the re-gauging induced 
by the transpositions (12), (13), (14) from the initial spanning tree to 
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Pl2 



Figure 3. Calculation of the action of (12) on T 3 

I \ /0 



the pushed forward one are 

/o i \ /o 



1 



A 



\ 



Pl3 



B*J 



A* 
1 



V 



,Pl4 



B 

c* 



V 



The calculation for p 12 can be read off from Figure 3. For this we read 
off the matrix $ from the re-gauging parameter and the matrix M a is 
given by the permutation we are considering. The other calculations 
are similar. All other transpositions, viz. those not involving 1, simply 
yield permutation matrices as there is no re-gauging involved. It is 
convenient to also have the following matrices as a reference: 



P(12)(34) 

and finally 



/0 1 
1 

V 



A 
B* 



P(14)(23) 



V 



B 



1\ 



c* 



P(123) 








i \ 


1 











A 





V 




B*j 



4.1.3. The point (0,0,0). At (0,0,0), the matrices pi2,P23,P34 give the 
usual representation of S 4 on C 4 . As is well known this representation 
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decomposes into the trivial representation and an irreducible 3-dim 
representation. This means that there is an at least 3-fold degenerate 
Eigenvalue A. Since the trace of H is zero, we also know that the 
Eigenvalues satisfy \x = — 3A. Plugging in (1, 1, 1, 1), which spans the 
trivial representation, we see that p = 3 and A = — 1. 

4.1.4. The point (tt, 7r,7r). In this case, the matrices pi2,P23,P34 only 
give a projective representation. As one can check that P12P23P12 = 
— pi 3 while P23P12P23 = P13 for instance. Define the 1-cocycle A by 
A(cr) = (—1) if 1 appears in a cycle of length > 1 and 1 else. So 
that A((12)) = A((13)) = A((123)) = -1 while A((23)) = A((24)) = 
A ((234)) = I. Then one calculates that p := po A has a trivial cocycle c 
and thus p is isomorphic to a true linear representation of § 4 . Checking 
the characters, one sees again that in this case the irreducible compo- 
nents of p, which also commute with H is again the one-dimensional 
trivial representation and the 3-dimensional standard representation. 
The trivial representation is spanned by (—1, 1, 1, 1). The Eigenvalues 
are then readily computed to be 1 with multiplicity 3 and —3 with 
multiplicity 1. 

Remark 4.2. We would like to remark that the choice of A amounts 
to choosing a different gauge for the root vertex, namely —1 instead of 
1. 

Remark 4.3. Notice that already in this case, even though there is no 
projective extension, our enhanced gauge-group is necessary. Without 
it there would only be an S3 action, those elements which involve no 
re-gauging. This smaller symmetry group is, however, not powerful 
enough to force the triple degeneracy, as §3 has no irreducible 3-dim 
representation. 

4.1.5. The point (tt/2, tt/2, tt/2) and{ -tt/2, -tt/2, -tt/2) . These points 
are similar to each other. We will treat the first one in detail. Again, 
we have only a projective representation of A4 aka. the tetrahedral 
group T. Namely, P(i 2 )(23)P(i3)(24) = -ip(i4)(23)- Again we can scale 
by a 1-cocycle A. This time \(id) = 1, \((ij)(kl)) = i, \{ijk) = 1 
if 1 ^ {i,j,k}, and X((ijk)) — i if 1 G {i,j,k}. The resulting repre- 
sentation p = p o A is then still a projective representation, but is it a 
representation of the unique non-trivial Z/2Z extension of A4, which 
goes by the names 2T, 2At, SL(2, 3) or the binary tetrahedral group. 
This group is well known. It is presented by generators s and t with 
the relations s 3 = t 3 = (st) 2 . In SX(2, 3) (that is the special linear 
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group of 2 x 2 matrices over the field with three elements F3), one can 

choose :;)aod ( =(:; 

For 2A 4 using a set theoretic section A of the extension sequence 
(13) 1 ^Z/2Z » 2A 4 ~^ A 4 > 1 

A 

and generator for Z/2Z, we can pick s = 2(123), i = z(234) 

as generators. Now we can check the character table, Table 1, and 
find that the representation p over the complex numbers decomposes 
as the sum of two irreducible two-dimensional representations X5 ®Xq- 
In fact, these are the two representations into which the unique real 
irreducible 4-dimensional representation of complex type splits over C. 
The explicit computation for the representation 

/ -i \ 



-% 
10 



(14) ~p{t) = -A((234))p (234 ) = - 



p(s) = -A((123))p (123) 

-Ij 

/l fl\ 
1 
10 
\0 1 0/ 

is as follows. Suppose the p = ©J =1 aiPi, where pi is the irrep with 
character Xi- Now tr(id) = 4, tr(— 1) = —4 , using the character table 
this implies that the coefficients a\ = a 2 = a.3 = a 7 = and furthermore 
(*) 04 + 05 + a6 = 2. We furthermore have that tr(s) = — 1 so that 
a 4 +ua 5 + w 2 a 6 = —1 which together with (*) implies that a 4 = 0, a 5 = 
a 6 = 1. This fixes the decomposition into irreps. As a double check 
one can verify that the rest of the equations are also satisfied. 

So indeed we find that (tt/2, tt/2, tt/2) is a point with two Eigenvalues 
with degeneracy 2. It is not hard to find (e.g. using the results of §4.1.6) 
that these Eigenvalues are ±v^3- 

The analysis of the complex conjugate point (—tt/2, —tt/2, —tc/2) is 
analogous. 

We would briefly like to connect these results to [3]. There it was 
shown that these four points are the only singular points in the spec- 
trum and that the two double crossing points are Dirac points. 

4.1.6. Super-Selection Rules and Spectrum Along the Diagonal. To il- 
lustrate the power of the super-selection rules we consider the action of 
the cyclic group generated by (234). One can easily see that the fixed 
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Representative 
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S 2 


T 


S 
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1 


1 


6 


4 


4 


4 


4 
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1 


2 


3 


3 


4 


4 


6 


Xi 


1 


1 


1 


1 


1 


1 


1 


X2 


1 


1 


1 


Ul 




a; 2 




X3 


1 


1 


1 


to 2 


uj 


uj 


uj 2 


X4 


2 


-2 





-1 


-1 


1 


1 


X5 


2 


-2 





— UJ 






u; 


Xs 


2 


-2 





-UJ 2 


— 




a; 2 


Xi 


3 


3 


-1 















Table 1. Character table of 2 • At [12], where u; = e s\ 



a, 6, c Group Iso class of type Dim of Eigenvalues A 







of extension 




Irreps 






(0,0,0) 


§4 


§4 


trivial 


1,3 


A = 
A = 


— 1 three times 
3 once 


(7T,7T,7r) 


§4 


§4 


trivializable 
cocycle 


1,3 


A = 
A = 


1 three times 
—3 once 


/ 7T 7T 7T \ 

2' 2/ 

/ 3-7T 3-7T 3-7T \ 

V 2 ' 2 ' 2 / 




2A A 


isomorphic 


2,2 


A = 


±\/3 twice each 






extension 









Table 2. Possible choices of parameters (a, b, c) lead- 
ing to non-Abelian enhanced symmetry groups and de- 
generate eigenvalues of H 



point set in the T 3 is the diagonal t = (a, a, a). The matrices \P actu- 
ally give a bona fide representation of C 4 . This is the representation of 
C 3 given by cyclicly permuting the last three factors of C. The action 
decomposes into irreps as follows: C 4 = triv@triv@uj@(jJ- Where oj is 
the 1-dimensional representation given by p((123)) — oj — exp(27ri/3). 
The two trivial representations are spanned by vi = (1,0,0,0) and 
v 2 = (0, 1, 1, 1) while the representation oj is spanned by w = (0, 1,uj,uj) 
and oj by id. 

Although we cannot extract information about the degeneracies from 
this it helps greatly in determining the Eigenvalues, since there are 
two irreps with multiplicity one each giving a unique 1-dimensional 
Eigenspace for the Hamiltonian. Hence we immediately get two Eigen- 
values. Plugging w and w into H(t) one reads off 



(15) Ai = uiexp{ia) + uj exp(— ia) A 2 = o;exp(m) + cuexp(ia) 
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Figure 4. Spectrum of H along the diagonal in T 3 



The sum of the two trivial representations gives a 2-dimensional 
isotypical component. Therefore, we have to diagonalize H inside this 
Eigenspace. It is interesting to note that at the special points it is 
exactly this flexibility that is needed in order to allow for crossings. 

To determine the two remaining eigenvalues A3 and A4, we apply H 
to v— xv\ + yv 2 = (x,y,y,y). The eigenvalue equation H v = A v 
leads to the equations 3y = \x and x + y(exp(ia) + exp(-ia)) = \y, 
Fixing x = 3 this gives the quadratic equation A 2 — 2cos(a)A — 3 = 
which has the two solutions 



(16) 



A 



3,4 



cos(a) ± -\/ cos 2 (a) + 3 



This gives the spectrum along the diagonal which is given in Figure 
4. The calculation only involves the classical symmetries without re- 
gauging. 

In loc. cit. the authors also assert that numerically they only found 
singular points in the spectrum along the diagonal. The fact that 
the arising candidates for Dirac points are indeed such points and the 
analytic proof that indeed there are no other singular points in the 
spectrum is contained in [3]. 

4.2. The P case. There is nothing much to say here. There is only 
the root of the spanning tree which is unique. The S3 action permutes 
the edges and their weights. This yields the permutation action on the 
T 3 . There is no nontrivial cover and the Eigenvalues remain invariant. 

4.3. The D case. Here things again become interesting. Permuting 
the two vertices, we obtain eight fixed points if a,b,c G {1, —1}. The 

matrix for this transposition is ( !? \ J . This gives super-selection 



rules and we know that v 1 



1 0/ 

^1, 1) and t> 2 



-1,1) are Eigenvectors. 
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The Eigenvalues being 1 + a + b + c and — (l + a + 6 + c) at these eight 
points. 

We can also permute the edges with the §4 action. In this case the S3 
action leaving the spanning tree edge invariant acts as a permutation 
on (a, b, c). The relevant matrices however are just the identity matri- 
ces and the representation is trivial. The transposition (12), however, 
results in the action (a,b,c) >->■ (a,ab,ac) on T 3 , see Figure 3. So to 
be invariant we have a — 1, but this implies that is the identity 
matrix. Invariance for (13) and (14) and the three cycles containing 1 
are similar. But, if we look at invariance under the element (12) (34) 
we are lead to the equations 

a = a,b = ac,c = ab 

This has solutions a = 1, b = c, for these fixed points again we find only 
a trivial action. But for a = — 1, b = — c these give rise to the diagonal 
matrix diag(l, —1) and hence Eigenvectors t\ = (1,0) and e 2 = (0, 1), 
but looking at the Hamiltonian, these are only Eigenvectors if it is 
the zero matrix if (a, b, c) = 0. Indeed the conditions above imply 
l + a + b + c = 0. 

Similarly, we find a Z/2Z group for (13) (24) and (14) (23) yielding 
the symmetric equations b = —1, a = —c and c = — l,a = —b. These 
are exactly the three circles found in [6]. Going to bigger subgroups 
we only get something interesting if the stabilizer group G t contains 
precisely two of the double transposition above. That is the Klein 
four group Z/2Z x Z/2Z. The invariants are precisely the intersection 
points of the three circles given by a = b = — 1 and c = 1 and its cyclic 
permutations. 

These are three lines of double degenerate Eigenvalue 0. These are 
not Dirac points since there is one free parameter and hence the fibers 
of the characteristic map of [3] are one dimensional which implies that 
the singular point is not isolated. 

4.4. The honeycomb case. This is very similar to the D story. The 
vertex interchange renders the fixed points a = ±l,6 = ±l which have 
Eigenvectors vi, t>2 as above and Eigenvalues 1 + a + b and — (1 + a + b) 
respectively. The irreps of the C3 action are triv © u. 

As far as the edge permutations are concerned the interesting one is 
the cyclic permutation (123) which yields the equations 

a = b,b = ba 

for fixed points. Hence a 3 = 1. We get non-trivial matrices at the 
points (ou,Cu) and (ui,u). At these points ei,e 2 are Eigenvectors with 
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Eigenvalue and H = 0, since 1 + a + b = 1 + uj + Co = 0. They are 
exactly the Dirac points of graphene. 

5. Conclusion 

By considering re-gaugings, we have found the symmetry groups 
fixing the degeneracies of the PDG and honeycomb families of graph 
Hamiltonians. The symmetries we used were those induced by the 
automorphisms of the underlying graphs. In our specific examples, all 
the graphs were highly symmetric, and hence had large automorphism 
group. Here we stress that our symmetries are extended symmetries 
and not just the classical ones. The most instructive and interesting 
case is the action of the binary tetrahedral group giving rise to the 
Dirac points in the Gyroid network. Note that as dimension-0 objects, 
the Dirac points for the Gyroid are codimension-3 defects in T 3 , rather 
than codimension-2 defects in T 2 , the case for the honeycomb lattice, 
which describes graphene. Nevertheless, one may expect that they too 
lead to special physical properties. 

There are several questions and research directions that tie into the 
present analysis. 

It would be interesting to find concrete examples of lifts of re- 
gaugings either in the non-commutative case or in the case of re- 
gaugings not induced by graph symmetries. One place where we intend 
to look for the former is in the non-commutative case of PDG and the 
honeycomb as we aim to probe the non-commutative/ commutative 
symmetry mentioned in Remark §3.7. 

We are furthermore interested in how these symmetries behave under 
deformations of the Hamiltonian and if they are topologically stable. A 
physically important type of deformations are those corresponding to 
periodic (in space) lattice distortions that describe crystals with lower 
spatial symmetry than those considered here. Such distortions may 
occur for instance during synthesis of the structure [1]. Codimension-3 
Dirac points, such as those of the Gyroid network, are especially in- 
teresting in this respect: they can be viewed as magnetic monopoles 
in the parameter space [13] and as such are topologically stable. This 
makes the physics associated with such points immune to periodic lat- 
tice distortions. 

Finally it seems that on the horizon there are connections between 
our theory and two other worlds. The first being quiver representations 
in general and the second being cluster algebras. The connection to 
the first is inherent in the subject matter, while the connection to the 
second needs some work. The point is that in our transformations, 
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we change several variables at a time. Nevertheless, the re-gauging 
groupoid can be viewed as a sort of mutation diagram. We plan to 
investigate these intriguing connections in the future. 
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